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"*^ ' Abstract. Using Monte-Carlo simulations on large lattices, we study the effects 

of changing the parameter u (the ratio of the adsorption and desorption rates) 

of the raise and peel model. This is a nonlocal stochastic model of a fluctuating 

interface. We show that for < n < 1 the system is massive, for u = 1 it is 

' ^ ■ massless and conformal invariant. For u > 1 the conformal invariance is broken. 

C ' The system is in a scale invariant phase with the dynamical critical exponent 

Q ' z and other critical exponents varying continuously with u. As far as we know 

O it is the first example of a system which shows such a behavior. Moreover in 

the broken phase, the critical exponents vary continuously with the parameter u. 
This stays true also for the critical exponent r which characterizes the probability 
v^ I • distribution function of avalanches (the critical exponent D staying unchanged) . 
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\^ ' The raise and peel model (RPM) of a fluctuating one-dimensional interface was 

^^ , introduced in Ref. [1]. In this model, the evolution of the interface (described in terms 

of RSOS paths), in the presence of a rarefied gas of tiles (tilted squares or particles), 

follows Markovian dynamics. When a tile hits a site of the interface it can be reflected, 

adsorbed with a probability Ua ("raising" the interface) or trigger a desorption process 

i-A . with a probability Ud, depending on the heights of the two neighbouring sites. In the 

f~| ' desorption process part of the top layer of the interface evaporates (the interface is 

O ■ "peeled"). The number of tiles which leave the interface in one desorption process, 

J-^ ' which may be macroscopic, define the size of an avalanche. 

The model depends on a single parameter u = Ua/ud or alternatively w = l/u. 
Lyj ■ The case u = 1 was discussed in detail in [1] in the continuum time limit. There are 

''^ , two reasons why this case has received so much attention: 

C^ ' a) The time evolution of the system is given by a Hamiltonian H which, in the 

finite-size scaling limit, has a known spectrum. The reason is that in this limit H 
describes a c = conformal field theory (CFT) (c is the central charge of the Virasoro 
algebra) and the spectrum is given by characters of this algebra [2]. Moreover this 
observation implies that the dynamical critical exponent is z = 1. 



C^ 



1. Introduction 
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b) The weights of the various RSOS paths in the probabihty distribution function 
(PDF) of the stationary state for a finite system are related to a counting problem 
of alternating sign matrices with different topologies (see [3] and references therein). 
Alternating sign matrices have miraculous combinatorial properties [4] and this has 
allowed to guess (conjecture) the expressions of several physical quantities in the 
stationary state for any number of lattice sites based on exact calculations for small 
lattices. 

In [1] very little work was done for u different of one except for a crude estimate of 
the phase diagram. The reason is simple: for u = 1 one has used the exact knowledge 
of the stationary state for lattice sizes up to 18 sites (it is hard, although possible, to 
go beyond this number) and made (unproven) conjectures for any lattice sizes. For 
u j^ 1, one has only diagonalized numerically (up to 14 sites) the Hamiltonian, that 
gives the time evolution of the system, in order to estimate in which domain of u the 
system is massless (scale invariant) or massive (finite correlation lenghts). The results 
were inconclusive. 

In this paper we present an extensive study of the model using Monte Carlo 
simulations. What was our motivation? 

a) For u = w — 1, the RPM is, to our knowledge, the first example of a stochastic 
process in which conformal invariance, a time-space symmetry, plays a role and one is 
interested to see some consequences. Let us stress that this is different of the one time- 
two space abelian sand-pile model where conformal invariance is seen in the stationary 
state [5] We also wanted to check the conjectures made in [1] for large lattice sizes and 
to study other properties of the model. 

b) When u — 1/w ^ 1 several scenarios are possible. The system may stay 
conformal invariant, may get massive or stay scale invariant without being conformal 
invariant. According to the standard lore (see Refs. [6] and [7] for a more precise 
definition of the conditions) provided that the interactions are local, scale invariance 
implies conformal invariance. Although for w = w = 1, the RPM can be mapped onto 
the spin zero sector of an Uq{sl{2)) invariant XXZ quantum chain [1] in which the 
interaction is local, as far as we can see, for u 7^ 1 a mapping onto a system with local 
interactions is not possible. Therefore, keeping in mind that the desorption process 
is not local, for u 7^ 1 we still have the interesting possibility that the system can be 
scale invariant without being conformal invariant. 

We will show that for m < 1 the system is massive while for m > 1 {w < 1) 
the system is scale invariant but not conformal invariant, with the critical exponent 
z [8] and other critical exponents varying continuously with w. At w = the system 
becomes massive again. We think that it is for the first time, that such a phenomenon 
is seen. 

c) As noticed already in Ref. [1] and confirmed in the present paper, for u — I 
the probability distribution function (PDF) of avalanches triggered by tiles hitting the 
interface in the stationary state, has a long tail, the critical exponents being D = \ 
and r = 3 (see Section 5 and Ref. [11] for definition of the exponents), which means 
that the model exhibits self-organized criticality (SOC). We will show how the PDF 
of avalanches behave if one varies w. 

The paper is organized as follows. In Section 2, following [1] we present the model 
and define the quantities which we found useful to characterize the interface. Those 
are contact points, clusters, heights and the fraction of sites of the interface where 
desorption does not take place (FND), i.e. the density of local maxima or minima. 

In Section 3 we describe the properties of the stationary states for various values 
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Figure 1. A configuration of the interface with three contact points and two 
clusters for the lattice size L = 16. 



of u, and derive the phase diagram of the model (Fig.^. We give the variation of the 
number of clusters with the size of the system L. We show that for u < 1 the density 
of clusters is finite for large values of L. For u = 1 this density vanishes like L~3 . For 
w < 1 the density of clusters also vanishes as a power law with an exponent which 
varies with w. We also give the local density of contact points as a function of the 
size of the system and its finite-size scaling behavior. We discuss the size dependence 
of the average height. The variation of the FND with w and the system sizes is given 
in Section 5. 

In Section 4 we give the values of the dynamical critical exponent z for various 
values of w. We estimate z using the Family- Vicsek scaling for various observables. 
We show that z decreases from the value 1 for w = 1 towards at small values of 
w. For w = 1 we show how the knowledge of the finite-size spectrum (explained in 
appendix A) gives information on the Family- Vicsek scaling functions. We also give 
the two-contact points correlation function obtained using Monte-Carlo simulations 
and illustrate the predictions of conformal invariance. 

In Section 5 we study the PDF of avalanches for different values of w. The average 
size of the avalanches stays finite in the whole domain 1 > u; > and its value is given 
remarkably well by mean field. We show that the exponent D is equal to 1 for all 
values of w but that the exponent r [11] changes from the value 3 at u; = 1 to a value 
close to 2 for smaller values of w. Our conclusions are presented in Section 6. 

2. The raise and peel model 

We consider a one-dimensional lattice with L + 1 {L = 2n) sites. An interface is formed 
by attaching at each site non- negative integer heights hi (i = 0, 1, . . . , L) which obey 
the restricted solid-on-solid (RSOS) rules: 

h,+i -h,^ ±1, ho = hL= 0, h, > 0. (1) 

There are C„ = (2n)!/((n -I- l){n\)^) possible configurations of the interface. In 
Fig. n we show a configuration for n = 8 (L = 16). Alternatively, one can describe 
the interface using slope variables Si = (/li+i — /ii_i)/2, {i ~ 1, .... L — 1). 

The dynamics of the interface is described in a transparent way in the language 
of tiles (tilted squares) which cover the area between the interface and the substrate 
{h2t = 0, /i2j+i = 1, (J = 0, ...,n)) (see Fig.^J. 

We consider the interface separating a film of tiles deposited on the substrate 
from a rarefied gas of tiles. We are interested to find the evolution of the interface 
towards the stationary state and to study the properties of the interface in this state. 

The evolution of the system in discrete time (Monte-Carlo steps) is given by 
the following rules. With a probability Pi = 1/{L — 1) a tile from the gas hits site 
i, (i = 1, ...L — 1). Depending on the value of the slope Si at the site i, the following 
processes can occur: 
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Figure 2. A dcsorption event. The incoming tile at site 1 triggers, with 
probabihty u^, an avalanche of 5 tiles, which are shaded. All of the shaded 
tiles are removed in the desorption event. 



i) Si ~ and hi > ft-i_i. 

The tile hits a local peak and is reflected. 

ii) Si = and hi < hi^i. 

The tile hits a local inininium. With a probability Ua the tile is adsorbed 
{hi i-^ hi + 2) and with a probability 1 — Wa the tile is reflected. 

iu) Sj = 1. 

With probability Ud the tile is reflected after triggering the dcsorption of a layer 
of tiles from the segment (hj > hi — hi^b, j = i + 1, . . . , J + 6— 1), i.e. hj i—f hj —2 
for j = i + 1, ...,i + b — 1. This layer contains 6—1 tiles (this is always an odd 
number). With a probability 1 — Ud, the tile is reflected and no desorption takes 
place. For an example see Fig. El 

iv) Si = -1. 

With probability u^ the tile is reflected after triggering the dcsorption of a layer 
of tiles belonging to the segment {hj > hi = hi^h, j = i — b + 1, . . . ,i — 1), i.e. 
/ij 1-^ /ij — 2 for J = I — 6 + 1 , . . . , I — 1 . With a probability 1 — Ud the tile is reflected 
and no desorption takes place. 

In the Monte Carlo simulations we choose Ud — i and Ua — u for u — Ua/ud < 1, 
and Ud — w,Ua = 1 for u > 1. For convenience we will use either u or w as the free 
parameter of the problem. We use these notations since we consider below either very 
small values of w (or u) and our intuition can be helped by the fact that for u = or 
w; = the properties of the system are known. Namely, for both u — and w — the 
spectrum is massive [1], the stationary states being the substrate {u = 0) and a full 
triangle {w — 0). 

It is useful to define some quantities which characterize the surface. An obvious 
geometric observable is the set of sites j {j even) for which hj = (the sites and 
L always belong to this set). These sites are also called contact points. This set is 
important, for example, to study desorption. Desorption events are limited to the area 
between two contact points, defined as a cluster (see Fig.^. As is going to be shown 
in Section 3, the density of contact points gives a local observable for which one can 
define various correlation functions (see also Section^. 

We first define quantities in the stationary state. The average number of clusters 
k{L) and the density of clusters p{L) are defined by: 

L 

k{L) = {J2S,^,o), (2) 

p{L) = k{L)/L. (3) 
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Figure 3. Typical configurations in the stationary states for L = 128 sites and 
tiiree values of u. 



Notice that the maximum value of p{L) is ^ which is obtained when u — and one 
has only one configuration in the system - the substrate. The average height is defined 
as: 



h{L) 



L 



(4) 



The lowest value of h{L), corresponding to the substrate, is ^. It is useful to consider 
also the average height in the middle of the system: 

h.{L)^{h.). (5) 

A relevant quantity is the average of the fraction of the interface where desorption 
does not take place (FND) 



n{L) 



1 



L-l 






(l-|-^.l)>- 



(6) 



This quantity will be used in Section 5 since it allows to estimate the average number 
of tiles desorbed in avalanches. The large L behavior of any of these quantities, say 
h{L), will be denoted hoc'- 

h^ = lim hiL). (7) 

L — ^cjo 

In Section 4 we are going discuss average quantities which are not only L 
dependent but also time dependent. For these quantities we use, for example, the 
notation k{t,L) to denote the number of clusters. 



3. The stationary states of the RPM 

In order to get an intuition on the behavior of the model when u changes, in Fig.|31we 
show for L = 128 typical configurations in the stationary states for three values of the 
parameter u. One notices that for u = 0.4 there are many clusters, fewer for u = 1 



Conformal invariance and its breaking 



I 



+ 



soc 



i 



u 



w 







Figure 4. Phase diagram of the raise and peel model. For < m < 1 (phase I), 
the model is massive. At u = nc = 1 it is massless and conformal invariant. For 
1 > ui > it is scale invariant with varying critical exponents and exhibits SOC. 
For UI = the system is massive. 
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Figure 5. The density of clusters poo for various values of u < 1. 



and a single one for w = 0.4 (other, less probable, configurations do exhibit several 
clusters). We also notice that the average heights are quite small. The Monte Carlo 
simulations were done for various lattice sizes up to L = 65536 sites. 

Before presenting our results, in Fig. 0] we show what be believe is the phase 
diagram of the model. In the domain < u < 1 the system is massive undergoing a 
second order phase transition ai u — Uc — I. For I > w > the conformal symmetry 
seen at w = 1 is broken but the system stays scale invariant. Since in this domain 
one has avalanches we call this phase a SOC phase. Another phase transition occurs 
at w = where the system becomes massive again. In the present section and in the 
next two ones we will give arguments which favor our conclusion. 

We first consider the density of clusters p{L) in the domain < u < 1. Taking 
large values of L, one obtains the u dependence of poo shown in Fig. [3 In all this 
domain the density of clusters stays finite. The density of clusters decreases from its 
maximum possible value 1/2 which corresponds to the substrate at u = to the value 
zero for u = 1. This indicates that there is a phase transition at u = 1. From our 
Monte Carlo simulations it was hard to obtain the critical exponent which gives p{u) 
when u approaches the value one. 

The local density of contact points g{x, L) for u — 0.4 and different lattice sizes 
is shown in Fig.El^a). One observes that for small values of x, and large values of L, 
g{x, L) decreases exponentially from the value 1 at x = (not shown on Fig.JHl^a)) to 
a constant value in the bulk %. 

X This value is twice that seen in Fig.|5] since for even x one always have g(x) = 0. 
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Figure 6. (a) The local density of contact points g{x, L) at u = 0.4 for different 
lattice sizes L. (b) The scaling function C(cos(7ra;/L)) defined in 1101 for u = 1. 



If we keep in mind that for u — one can diagonalize the related Haniihonian 
of the stochastic process and show rigorously [1] that the system is massive, we can 
expect that the system is massive in the whole domain < u < 1. 

We turn now to the interesting case u = 1. In Ref. [14] a conjecture was made 
which gives the probability to have k clusters in a system of size L. For our purpose, 
we need only the behavior of the average number of clusters (see [1] Eq.(47)) for large 
system sizes, given by 



k{L) 



r(i/3)V3 I 



(8) 



which shows that the density of clusters vanishes algebraically for m = 1. We have 
checked this conjecture using Monte Carlo simulations up to a system of size L = 512. 
Since, as we show below, the density of clusters also vanishes algebraically in the entire 
domain w < 1, it is convenient to define a critical exponent a < 1 which gives the 
power law increase of the number of clusters with the system size, k{L) ^ L". For 
u = 1 we have obviously a — 2/3. 

We have considered, again for u = 1, the density of contact points g{x,L) for 
several lattice sizes and noticed that for the function 

G(x,L)=iig(x,L) = G(y), (9) 

where y — x/ L one has a data collapse. We make an ansatz about the functional 
behaviour of G{y) taking: 



,L 



g{x,L) = [-sin(7ra;/L)]-i/^C(cos(7ra;/L)). 



(10) 



The motivation of this ansatz will become clear in the next lines. In Fig. |^b) 
one shows the function C(cos(7ra:/i). It is practically a constant. If one assume that 
C is independent of nx/L, one can get its value by integrating g{x, L) (Eq. ((Tn|l ) and 
compare the result with the average number of clusters given by (O. One obtains: 



C 



^r(4) 



0.753149. 



67r6 6 

in very good agreement with the data (see Fig. Etb)). 



(11) 
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Figure 7. Average number of clusters fc(I/) as a function L for various values 
of w . The values of the exponent a (slope) is also shown. 



If the density of contact points corresponds to a local operator in conformal field 
theory, one expects [15], 



g(x, L) = C[— sin(7ra::/L)] 



X 



(12) 



where X = A + A is the scaling dimension of the local operator (A — A is the 
conformal spin). Comparing Eqs. (|10|I - H12() one concludes that the density of contact 
points corresponds to a local operator with scaling dimensions X — 1/3. This result 
is surprising. One can show that the profile of an operator with spin vanishes [15]]. 
This implies A = A = 1/6. This value can't be obtained from the possible values of 
A (see (|A.8|l ) of minimal models. It might be obtained from iy(^„) algebras [16, 17] 
but it is not clear why W-symmetry should play any role in our problem. Another 
possible explanation is that for a conformal theory with c — (this is also the case of 
percolation), under certain circunstances, the scaling dimension can be an arbitrary 
number [18, 19]. This possibility is a bit strange since the number one has to obtain 
is a neat 1/6. Finally, it is possible, that the proof that operators with spin have no 
profile functions does not apply to our problem. If this is the case, one can choose 
A = and A — 1/3 §. The value 1/3 is a perfectly acceptable one (see Appendix A). 
To conclude, we have not yet an explanation for the value X = 1/3 in Eq. (|12(l . 

We move away from u — 1 and consider the domain 1 > w. We observe that in 
Fig.Elthe number of clusters k(L) keep having an algebraic increase with an exponent 
a that decreases monotonically with w. In Table^we give estiinates for the exponent 
a for several values of w. We have also checked for several values of w that the scaling 
law 



G(x,L)=Li-"g(x,L) 



G(f) 



(13) 



holds. 

From this investigation one can conclude that for values of w > 0.05 and maybe 
smaller, the system stays scale invariant with a varying critical exponent. In the next 
sections we are going to show that this picture stays valid. 

§ For the appearance of operators with conformal spin in the presence of the quantum group 
symmetry mentioned in Appendix A, see Ref. [20]. 
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1 


0.85 


0.75 


0.4 


0.25 


0.1 


0.05 


0.025 


a 


0.67 


0.50 


0.40 


0.30 


0.24 


0.15 


0.06 


0.01 



Table 1. Estimates of the exponent a giving the increase of the number of 
clusters with the size of the system L. 
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Figure 8. Average height h{L) as a function of the lattice size L for various 
values of to > 0.1. 
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Figure 9. (a) Average height in the middle of the lattice hi (L) as a function of 

2 

the lattice size L for some values of w in the range 1 > ui > 0.025. (b) Average 
height in the middle of the lattice hi{L) as a function of the lattice size L for 

■2 

various small values oi w {w < 1/200). 



Although it seems obvious that one should look at the L dependence of the average 
value of the heights h{L) or the average value of the height at the middle of the lattice 
hi{L), it turns out that these quantities are harder to study since they vary slowly 
with the system size. In Fig.|Slwe show h{L) for various values of w > 0.1. One sees 
a slow (logarithmic ?) increase when L increases. A similar picture is obtained if one 
considers hi(L) (Fig. Efa)). 

A different picture emerges if one considers very small values of w as shown in 
Fig. Efb). One can observe two different phenomena. The values of hi{L) increase 
substantially as compared to the values observed at larger values of w, and more 
interestingly, it looks like they saturate for large values of L. Two scenarios are 
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Figure 10. Average height in the middle of the lattice hi(L) as a function of 
w for the lattice size L = 65536. The dashed curve is given by I14i . 

compatible with our last observation. There is a change of the physics which occurs 
at very small values of w illustrated by a change of the behavior of the L dependence 
of hi{L). This is the first scenario. If this would be the case, we would have observed 
a crossing of the increasing functions hi (L) obtained for larger values of w and the 
flat values of hi{L) observed for small values of w. This is not the case. What looks 
more plausible is the second scenario where for larger values oiw, hi (L) also saturates 
but this is seen only at larger values of L than those considered in our simulations. 

According to the second scenario, for all values oi w < 1 (maybe including w = I) 
and very large values of L, hi{L) saturates. In order to check if this scenario makes 
sense, taking L = 65536 we have done a fit to the curve which gives the maximum 
values of hi{L) (denoted by /irnax) for various values of w. One obtains: 

• "' (14) 



tmax 



w 



,0.62 ■ 



As seen in Fig. 1101 one obtains a fair fit to the data. Interestingly, if we use (|14|l for 
w = 1 we obtain /imax = 4.67, a value compatible with the results shown in Fig.|^a). 
Taking the results shown in Fig. Ob) at face value, implies that the phase SOC (see 
Fig. ^ extends at least to the value w = 10~^. On the other hand for w = the 
single configuration is the full triangle, /irnax has values of the size of the system 
and therefore the L-dependence of /imax shows no saturation. This suggests a non- 
analytical behaviour of various quantities describing the system ai w — Q. In order 
to illustrate this point, in Fig. ^2 we show typical configurations in the stationary 
state for a very small value of w {w = 0.00025) and different system sizes. One can 
see a clear change in the role of the boundaries if L varies from 1024 to 32768. For 
L = 1024 one is close to the full triangle configuration (expected for w = 0). However 
as one increases L one gets a plateau (in which the height is almost constant and L 
independent) and the boundaries play a less and less important role. 



4. Space-time phenomena 

In the last section, in which we have studied the properties of the stationary states, 
the data have confirmed the phase diagram shown in Fig. 0] In the present section 
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Figure 11. Typical configurations in tlie stationary state for w = 0.00025 and 
diflferent values of L. 



we are presenting new facts which suggests that our assumption is correct. We are 
going to use the Family- Vicsek [21] scahng in order to check if the system is in a 
scale invariant phase and if the answer is positive, determine the dynamical critical 
exponent z. If we consider a time-dependent average quantity a{t, L), being a{L) the 
average in the stationary state, than for large values of t and L the function A{t^ L) 
scales as: 



ait, L) . , t ^ 



(15) 



where z is the dynamical critical exponent. The scaling function depends on the initial 
conditions. In the data shown below, we have taken as initial condition the substrate 
{h{2i) = 0, h{2i -I- 1) = 1, (i = 0, . . . , n)) with probability one. We first consider the 
quantity K{t, L) which corresponds to the average number of clusters {a{t, L) in (|15|l 
is k{t,L) in this case). 

In Fig. ll2r a) we show the function K(t, L) for u — 0.4 for different lattice sizes. We 
notice that for large lattice sizes, K{t, L) is a function on t only and can be obviously 
fitted with a sum of exponential functions as expected in a massive phase. Fig. I12r b'l 
shows K{t, L) for m = 1 and as expected from conformal invariance, z — 1. Conformal 
invariance gives also information on the function K[t, L). One can use the knowledge of 
the finite-size scaling limit spectrum and obtain the functional dependence of K{t/L). 
To illustrate how it works (see Fig. I12f b'l'l. for large values of t/L one can do a fit to 
the data, using the following ansatz: 



K{t/L) ^ Cie-— + C2e-— , 
where the Xi are obtained from the spectrum of the Hamiltonian ljA.7|) , 



E, 



L 



, z = l,2. 



(16) 



(17) 



The constants Ci depend on the initial conditions. If one starts with the substrate 
Ci = 0.28 and Ca = 0.54. 
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Figure 12. (a) Function K{t,L), defined on 1151 . for u = 0.4 as a function 
of time for tlie lattice sizes L = 2048, 4096 and 8192. (b) Function K{t,L) for 
M = 1 as a function of t/L for various lattice sizes. Red curve is given by a fit 
to 1161 considering only the first exponential with Ai = 7r^\/3 and Ci = 0.18. 
In the inset we compare K(t,l) — C\e~'^^^'^ for L = 8192 with the short time 



exponential C2e 



-Aat/i 



(given in red line) with A2 = 4Ai and C2 = 0.54. 



w 


K 


H 


Hl/2 


1 


1.03 


0.95 


0.97 


0.9 


0.86 


0.74 


0.72 


0.7 


0.55 


0.52 


0.50 


0.4 


0.35 


0.34 


0.34 


0.25 


0.31 


0.30 


0.30 


0.025 


* 


0.07 


0.07 



Table 2. Estimates of the dynamical critical exponent 2 for various values of 

w. The estimates were obtained using 1151 with A replaced hy K, H and Hi. 

2 
The number represented by (*) is not reliable since in that region the number of 
clusters is very small. 



We now consider smaller values of w (larger values of u). We have done an 
extensive study of various functions A{t,L) taking k{t,L), h{t,L) and hi{t,L) in 
H15|l . The estimated values for the critical exponent z are shown in Table 2. For very 
small values of w, the estimates coming from K{t^ L) are poor because the number of 
clusters is small and hence large errors. 

Inspection of TableElshows a smooth drop of the values of z from 2; = 1 for w = 1 
to z = for w = 0. The latter values is in agreement with a direct calculation of the 
mass gap at w = [1] which gives a finite gap in the thermodynamic limit. To show 
how the estimates given in Table-work, we demonstrate in Figs.^Ja) and ll3f b) the 
scaling functions K{t, L) for w — 0.4 and Hi (t, L) for w — 0.0025, respectively. 

We have not studied extensively the two-contact points space-time correlation 
functions except for the special case u — 1. This case is interesting since from 
conformal invariance the two-contact points space-time correlation function exhibits, 
far from the boundaries, the following scaling form in the bulk (see (O): 



C{R,t,L) = L-iG{fi), 



(18) 
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Figure 13. (a) K{t,L) for w = 0.4 as a function of t/L'S for the lattice sizes 
L = 1024, 2048, 4096 and 8196. (b) Hi(t,L) = hi (t, L)//ii(L)-l for «; = 0.0025 
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as a function of t for various lattice sizes. 
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Figure 14. Tlie 2-contact point correlation function C'{R,t,L) times La as a 
function of the scaled distance R/L and the scaled time ^4 (vg = §\/3 is the 
sound velocity) for various lattice sizes. 



where 



M 



L ' 



(19) 



In (|18|l and H19() R is the distance between the contact points, t is the time difference 
and Vs = |\/3 is the sound velocity (see Appendix A). The factor L~9 in H18|l is 
obtained from 0. In Fig. El we show the scahng function G{fi) ||. To avoid boundary 
effects the data was taken only in the center segment of size L/2. 

The functional dependence (|18|l which is special for the conformal invariant case is 
nicely seen in Fig. 1141 Our data are not good enough in order to extract from the small 
fj, behavior of G{^), the two-contact point correlation function in the thermodynamic 
limit. 



We have not used the connected correlation function in 1181 . 



Conformal invariance and its breaking 14 

5. Avalanches in the RPM 

An interesting property of the RPM is the occurrence of avalanches in the desorption 
processes. For the case u = 1, based on exact results on small lattices, the existence 
of avalanches was already suggested in [1]. In this section, using Monte Carlo 
simulations on large lattices, we present the properties of avalanches in the whole 
domain < w < 1. 

In the stationary state, once a tile from the rarefied gas hits the interface, it can 
be reflected, adsorbed or can trigger a nonlocal desorption process in which many tiles 
may leave the interface, this defines an avalanche. In the latter case the size of the 
avalanche is given by the number of tiles T that are released in the process. This 
number is always odd, therefore it is convenient to write: T = 2v — 1 (w = 1, 2, . . .). 
For w < 1 only a finite number of tiles are removed, since the density of clusters is 
finite. For w < 1, the cluster density vanishes, and therefore macroscopic number of 
tiles may be desorbed, hence macroscopic avalanches. 

We denote by S{v, L) the PDF which gives the probability for an avalanche of 
size V for a system of size L. In models of self-organized criticality (SOC) [5,11,22-24], 
for large values of v and L, one may expect the PDF to exhibit a simple scaling law, 

S{v,L)^v-^F{^), (20) 

characterized by two exponents r and D. 

It is convenient to consider moments of the PDF and their scaling properties: 

< «™ >L= J2 ""^("' ^) = r™-^"'"^ (21) 

v=l 

where 

/ N f 0, TO < T - 1 , , 

ct(to) = ^ (22) 

^ ^ [ D{m + 1 - t), to>t-1 ^ ^ 

In order to get estimates for D and r from the Monte Carlo simulations, it is 
convenient to take couples of system sizes L and L' and consider the quantities: 

aL,L' im) - In(f^)/ In(-) = mDL,L' + Bl,l' (23) 

(w™)l' L' 

and 

0'L,L'(wl,l') =0, Ti^L/ = TOL_i/ + 1. (24) 

Using the results of Monte Carlo simulations for couple of sizes (in practice we 
have chosen L' = 2L) one can obtain estimates for the two critical exponents D and 
T from aL,2L and tl^2L- In Fig.^|we illustrate the to, dependence of aL,2L obtained 
from Monte Carlo simulations, for several values of L in the u = 1 case. 

From a cursory look at this figure one can see that for u = 1, r sa 3 and £> w 1 
confirming the results obtained in [1] using much smaller lattices (for a more precise 
estimate, see Table EJ. The result D = 1 is to be expected since in a conformal 
invariant theory (this is the case if u = 1) one has no other length than the size of 
the system L. We found no explanation for the value three of r although this number 
shows up in many aspects of the model (see for example © and the Appendix). In 
Fig- El we show the scaling function F{v/L) defined in (|20(l . One observe a nice data 
collapse. 

In Table ini we give the estimates for the exponent r and D for various values of w. 
The estimates presented in this Table are obtained using only one pair of sizes (4096 
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Figure 15. The m dependence of c7l,2L fof w = 1. The fitted straight line is 
also shown. 



L = 1024 L=2048 




Figure 16. The scaling function F{v/ L) for ui = 1. The data are for lattice 
sizes L = 1024 and L = 2048. 



and 8192) but we made sure that the results are rehable by studying the changes of 
the values of the estimates using smaller lattices. 

For small values of w, the determination of r has to be done with care. This is 
due to surface effects (there are avalanches occurring at both ends of the system which 
have a different behaviour than avalanches occurring in the bulk). The phenomenon 
is illustrated in Fig. I17r a) where the PDF S{v, L) is shown for w = 0.00025 and 
L = 32768. One notices that S{v/L) has two peaks: one for v — \ and another one for 
very large values of v. Large avalanches occur when one tile hits the boundary of the 
interface. On the other hand, one expects that for large values of L, the probability 
for a tile to hit a boundary get smaller and smaller as the system gets larger (see 
Fig. Ill|l . therefore the probability to have large avalanches should decrease with the 
size of the system. This phenomenon is nicely illustrated in Fig. llTf b). where one can 
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Figure 17. (a) The PDF for avalanches for w = 0.00025 and L = 32768. (b) 
The PDF of avalanches for w = 0.00025 and different sizes. 
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Figure 18. (a) Estimates ri^^^L for w = 0.00025 using different values of L from 
L = 2048 to L = 16384 (b) The m dependence of crL^2L for w = 0.00025. 



see that the probabiHty to have an avalanche at the boundary decreases with the size 
of the system. One expects that the determination of r using (|23|l and (|24|l will give a 
smaller value for t for small lattices (large avalanches at the boundaries) and a larger 
value when compared with the value obtained by using large lattices (the importance 
of the boundary effects decreases). That is exactly what happens as can be seen in 
Fig. Iisr a). One sees in this figure that the estimates for r converge to a value close 
to 2. 

In order to illustrate the convergence of the estimates in Fig. llHr b) we present 
for a very small value of w the functions (Tl.2L for several values of L. Although the 
convergence is poorer when compared with these quantities obtained for higher values 
of w (see Fig. ll5|) . one can still get reliable results. 

The value of the exponent D stays unchanged and equals to one in the whole 
range w < 1. This can be understood in the following way [25] the tiles which are 
desorbed and create the avalanche belong always to a one-dimensional layer. 

We turn now to the exponent r. As can be seen in Tabled this exponent varies 
from the value 3 (u = 1) to a value close to 2 for w close to zero. In the case u = 1, 
the average size of the avalanches stays finite when L gets large but the dispersion 
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w 


1 


0.95 


0.90 


0.80 


0.60 


0.45 


0.30 


0.05 


0.02 


0.005 


Dl,L' 


1.004 


1.089 


1.087 


1.066 


1.040 


1.026 


1.017 


1.002 


1.002 


1.006 


tl.u 


3.00 


2.77 


2.63 


2.47 


2.32 


2.25 


2.18 


2.07 


2.046 


2.00 



Table 3. Estimates of the critical exponents D and r obtained by using lattices 
sizes L = 4096 and L' = 8196. 
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Figure 19. The functions F{v/L), defined on l2Ui . at w = 0.00025 for some 
values of L. 



diverges (the function F{v/L) shown in Fig.^ldoes not vanish at the origin). Does an 
exponent r = 2 mean that the average size of the avalanches diverges logarithmicaUy 
with L ? Not necessarily, it depends if the scaling function F{v/L) vanishes or not at 
the origin. In Fig. ^]we show for w = 0.00025, that F{v/L) does not vanish at the 
origin. Since later in this section we are going to show that (v) stays finite for any 
w > 0, we conclude that r gets very close to the value 2 but never reaches it. 

We would like to comment now on the properties of the avalanches seen at the 
boundaries and illustrated in Fig. 17. For large values of v where the boundary effects 
occur, the PDF S{v, L) presents two striking features: the probability to have an large 
avalanche decreases with L and S{v,L) increases with v. A close inspection of the local 
heights distribution for different values of L seen in Fig. E] allows us to investigate 
both features. Firstly, one sees that for large values of L the boundary becomes 
less and less important and therefore the probability for a tile to hit the boundary 
and produce a large avalanche decreases. On the other hand, at the boundary, the 
probability for a hit by a tile is proportional to \dh(x, L)/dx\ (keep in mind that a 
desorption process takes place when the local slope s{i) takes the values ±1). From 
the data shown in Fig.llllone can see that for x close to the boundaries, d'^h{x, L)/dx^ 
is negative. This implies that the probability to hit the interface with a tile at the two 
ends of the system is larger than further into the system. 

An important consequence of this study is that avalanches exist at least in 
the interval 1 > w > 0.00025 and therefore the system is scale invariant in this 
interval. (The data shown in Tables ^ and |21 are limited to the interval 1 > w > 0.025). 
It is probably safe to assume that the SOC phase covers the domain 1 > w > 0. (We 
have obviously checked that for u < 1, where the system is massive, for large values 
of V the PDF has an exponential decrease and not an algebraic one) . 

We turn now our attention to the average size of the avalanches. A simple mean- 
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field argument allows to compute it using informations about the stationary state 
only. 

We first consider, in the stationary state, the average fraction of interface where 
desorption does not take place n{L) (FND) defined in IjHJi. This quantity varies from 
the value 1 for the substrate (the single configuration for m = 0) to the value 1/(L— 1) 
for the full triangle (the single configuration for w = 0). 

For M = 1 there is a conjecture [26] for the values of n(L) 

This conjecture was checked using our Monte Carlo simulations for various lattice 
sizes. 

If one knows n(L), the average probability Pa{L) to have an adsorption process, 
the average probability Pr{L) to have a reflection process and the average probability 
Pd{L) to have a desorption process can be easily obtained: 

P^(i) ^ !l(^ _ _i_, (26) 

^ ^ 2 2(L-1)' ^ ' 

Pr{L) = ^{2w-l) + l-w + ,^^^—^y (27) 

Pa{,L)+Pr{L) + Pd{L)^l, (28) 

where, like in our Monte Carlo simulations, we choose Ua — ^ and Ud = w = \/u. 

Taking into account that in the stationary state the average number of adsorbed 
tiles is equal to the average number of desorbed tiles, one obtains: 

(T)i = 2(^;>i-l^^. (29) 

Notice that this is a mean-field calculation since we have first computed the average 
probability to have a desorption process and multiplied it with the average number of 
tiles which are desorbed. In the large L limit 126|) - (|29|l gives: 

{T)oo « ,,. "°° , . (30) 

2(1 - noojui 

In Fig. [201 one shows the w dependence of noo as obtained from the extrapolated 
results of our Monte Carlo simulations. One notices a discontinuous behavior of rioo 
around w = Q. At w = one has Uoo = 0, while limu,_^o+ '^oo = 0.5. Using H30|l one 
concludes, in agreement with previous observations that for any finite value w, the 
average size of the avalanches is finite. 

In Table 21 we compare in two cases the mean- field predictions for (r)^ given in 
H3U|I with the values measured in Monte-Carlo simulations. The agreement is excellent. 
In the same Table the values for n{L) for several values of L are also given. One notices 
that for u ~ \ and L = 1024 one already gets within five digits the asymptotic value 
3/4. This is not the case for w = 0.00025, where the finite-size effects are larger. Our 
numerical analysis indicate that the asymptotic behavior of n{L) can be described by 

n(L) = noo(l-A(w)/L + ---), (31) 

where ^(1) = 1/6 and A{'w) diverges if w; ^ (data not shown). 

We conclude this section with the following observation, which is necessary for 
the consistency of the present model for avalanches. If one looks at the expression 
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Figure 20. The w dependence of the fraction of the interface where adsorption 
does not take place rioo obtained from the extrapolated results of Monte Carlo 
simulations. 



w 


L 


UL 


{T)l 


MF 


{T)l /MF 


1 


1024 


0.749 


1.497 


1.497 


0.999 


1 


2048 


0.749 


1.498 


1.501 


0.998 


1 


4096 


0.749 


1.499 


1.500 


0.998 


1 


8192 


0.749 


1.499 


1.499 


0.999 


0.00025 


1024 


0.199 


496.999 


496.119 


1.001 


0.00025 


2048 


0.313 


914.346 


913.820 


1.001 


0.00025 


4096 


0.407 


1374.713 


1373.856 


1.001 


0.00025 


8192 


0.455 


1674.877 


1674.550 


1.000 



Table 4. Values of rioo and the average number of tiles obtained from Monte 
Carlo simulations (r)^ as compared with the compared with the mean- field results 
(MF). The values are given for several lattice sizes L and two values of ui, ui = 1 
and w = 0.00025. 



H27|l of the average (over configurations) probability that a tile hitting the interface 
is reflected, one sees that for small values of w, the probability approaches the value 
one. This implies that to obtain large avalanches (small r) one has to wait a long 
time. This is to be expected in models of SOC. 



6. Discussions 



The raise and peel model describes a stochastic process of a one-dimensional interface. 
RSOS paths give the configuration space. It is well known that if a stochastic model 
has as stationary state in which the RSOS paths have the same probability, the average 
height {h{L)) increases with the system size L like vL (Brownian motion). In the 
one-parameter (u or it; = l/u) RPM, in the stationary states, the various RSOS paths 
get weights and henceforth the physics is different. It turns out that for all finite 
values oi u {w y^ 0) in the large L limit, {h{L)) stays finite. The convergence to the 
large L limit of {h{L)) is very slow. The w = value is singular: in the stationary 
state one has, with probability one the full triangle and therefore h — L/2. The phase 
diagram is shown in Fig. ^ For m < 1, one is in the massive phase (finite correlation 
lengths and non- vanishing densities of clusters). For 1 > w > the density of clusters 
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vanishes algebraically with an exponent which depends on w (Table^. The dynamical 
exponent z also depends on u. It varies between one and zero (see TableEJ- We called 
this phase the SOC phase since in the desorption process one has avalanches. The 
PDF of the number of tiles desorbed has long tails. Although the average number 
of tiles involved in the avalanches stays finite for large L, higher moments diverge. 
The r exponent changes from the value 3 (w = 1) to the value 2 for w close to zero 
(Table Oil. The exponent D is robust and keeps the value D = 1. A fascinating aspect 
of the model is that, contrary to our intuition, for any small value w > 0, if L is 
large enough, the system separates into a bulk part and a "surface" part related to 
the condition that at the two boundaries the heights have to take the value zero (see 

Fig. nil). 

The value w = 1 is special. In the continuum time limit the Hamiltonian which 
gives the Markovian time evolution of the stochastic process, has the same spectrum 
as the XXZ quantum spin 1/2 chain Hamiltonian with a given Z-anisotropy and fixed 
boundary terms (see Appendix A). This Hamiltonian is integrable and the finite-size 
scaling spectrum is known. The main observation is that for u = 1 one has a space- 
time symmetry (conformal invariance). There are several important consequences 
concerning the properties of the system. The one-point function which gives the local 
density of contact points has the scaling form given by H12|l . The two-contact points 
space-time correlation function has a scaling form given by (|18|l . The dynamical 
critical exponent z can be computed analytically and has the value 1. There are 
other consequences of conformal invariance. The Family- Viczek scaling function which 
gives the time and size dependence of average quantities can also be obtained (see 
Fig. ll2f b)). The value D = 1 is a consequence of the fact that in a conformal theory, 
in the large L limit the single scale is L itself. The value r = 3 remains a puzzle 
although the number 3 appears all over for u — 1 (see Appendix A) . 

The reader not familiar with conformal invariance might skip the end of this 
section. We have shown that the density of contact points gives a local operator of 
a c = (c is the central charge of the Virasoro algebra) conformal field theory. As 
discussed in Section 3 the value of the scaling dimension of this operator has not a 
simple explanation. Other local operators can be defined in a similar way. For example 
one can look at the local density of points where one has a deep (slope zero) at the 
height 1 (or any fixed height). 

One interesting aspect of the RPM is that one has a simple example of an unusual 
kind of conformal symmetry breaking. For u < 1 the system becomes massive. This 
is not a surprise. One can understand why this happens in the following way: one 
can take the desorption rate Ud = 1 (see Section|21l and Ua = u. This implies that the 
perturbation is local (the adsorption is local), li w < 1 the situation is quite different 
and here comes the surprise. One has scale invariance without having conformal 
invariance. The appearance of this novel phenomenon has the following explanation. 
For w < 1 one takes Ua = 1 and Ud — w. Desorption is a nonlocal process and 
therefore the conformal theory is perturbed nonlocally and there is no a priori reason 
why conformal invariance cannot be broken while scale invariance is not. 
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Appendix A. The finite-size scaling limit of the spectrum of the 
Hamiltonian for u = 1. Results from conformal field theory. 

As shown in Ref. [1], in continuum time, the time evolution of the RPM for u = 1, is 
given by the Hamiltonian: 

L-l 

H=J2{l-e,), (A.l) 

1=1 

where e^ acts on the configurations (RSOS paths) described in Section 1, taking one 
path into another or leaving it unchanged. Using a similarity transformation from the 
configuration basis to a spin basis, the e^'s have a representation in terms of Pauli 
matrices: 

e^ = lwf<T^+i + ^f ^r+1 - l<y^<yUi + 1 + *^« - <+i)]- (A.2) 

If we introduce the expressions ljA.2|) in (jA.l|l . one obtains a well known 
Hamiltonian of an integrable quantum chain acting into a 2^ dimensional vector space. 
This chain has Uq{sl{2)) symmetry with q = exp(«7r/3). To obtain the action of H in 
the vector space of the RSOS paths one has to consider only the subspace of Uq{sl{2)) 
scalars (spin s = 0) which has the dimension C„ = (2?T,)!/(n+l)(n!)^ with L = 2n. The 
spectrum of the Hamiltonian IjA.ip is non-negative. Let us denote by Er (r = 0, 1, . . .) 
the energy levels in increasing order: 

Eo = 0< Ei< E2< ■■■. (A.3) 

The partition function giving the finite-size limit of the spectrum of the Hamiltonian 
of the quantum chain, is defined as follows: 

Ziq) = Km ZL{q) = lim ^ q^ , (A.4) 

where Vg — 3-\/3/2. One can show [13] that Z{q) has the expression 

OQ 

Z(g) = 5](2s + l)e«((?), (A.5) 

s=0 

where 

oo 

6(g)-<7'^^^(l-g'^+^)n(i-'?")~' (A-6) 

n=l 

is the partition function in the spin s sector. 

In Section 4 we look at the time evolution of the stochastic model. The finite 
size-scaling limit is given by ()A.6p with s = 0. This implies that in the large L limit 
we have 

LEo = 0, LEi = 7r3\/3, LE2 = t^i^Vs, (A.7) 
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Up to now we have described the time evolution of the system. If one considers 
space-time correlation functions, the various critical exponents are expected to be 
obtained [27, 28] from the scaling dimensions of the primary fields of minimal models 
[29] with a central charge c=0 [1]: 

where p and q are positive integers. Notice that the value A = 1/3 appears in both 
the finite-size spectrum of H (see (|A.6|I with s = 1) and in l|A.8|l . 
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